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Abstract 



The construction of the cylinder at spatial infinity for stationary spacetimes is considered. 
Using a specific conformal gauge and frame, it is shown that the tensorial fields associated to 
^-l the conformal Einstein field equations admit expansions in a neighbourhood of the cylinder 

at spatial infinity which are analytic with respect to some suitable time, radial and angular 
coordinates. It is then shown that the essentials of the construction are independent of the 
choice of conformal gauge. As a consequence, one finds that the construction of the cylinder 
at spatial infinity and the regular finite initial value problem for stationary initial data sets 
are, in a precise sense, as regular as they could be. 



> 

^ ' 1 Introduction 

O 



The present article discusses a certain conformal extension for vacuum stationary solutions for 
the Einstein field equations — the so-called cylinder at spatial infinity. This construction provides 
detailed information about the structure of this class of spacetimes in the region where spatial 
infinity touches null infinity. The analysis presented here is key for the construction of asymptot- 
ically simple spacetimes from initial value problems on Cauchy hypersurfaces. In order to better 
understand the context of our analysis, we present a brief overview of the ideas and problems 
involved. 

1.1 Asymptotically simple spacetimes 

Penrose's notion of asymptotic simplicity — see e.g. [351 136] was introduced with the objective 
of providing a framework for the discussion of isolated systems in General Relativity. The pro- 
gramme behind this idea is usually know as Penrose's Proposal — see e.g. [23 HI] ■ A vacuum 
spacetime {M,g^ v ) is said to be asymptotically simple if there exists a smooth, oriented, time- 
oriented, causal spacetime (A4,g M „) and a smooth function S (the conformal factor) on M. such 
that: 

(i) M. is a manifold with boundary ^ = dM.\ 

(ii) e>0onM\/,ande = 0,de^0 on J; 

(iii) there exists an embedding <f> : M \ J — >• M such that g^ v — 2 ((f> _1 )*5 /J!y ; 
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(iv) each null geodesic of (M,g^) acquires two distinct endpoints on ,P . 

In this definition, as well as in the rest of the article the word smooth will be used as synonym of 
C°°. In order to simplify the notation we will write g^ v = Q 2 g^ instead of g^ v = 2 ($ _1 )*5 All/ . 
The point (iv) in the definition excludes black hole spacetimes — however, the discussion in this 
article will be local to the conformal boundary ^ ', and hence (iv) will not be of relevance in our 
considerations. 

The first natural examples of asymptotically simple spacetimes are solutions to the Einstein 
field equations which are static or, more generally, stationary near the conformal boundary. That 
this is the case is a consequence of the structural properties of the static and stationary field 
equations — see e.g. [T5]. We shall elaborate further on this in the coming paragraphs. Now, 
in order to be of real physical value, the notion of asymptotic simplicity should also include 
spacetimes which are not static or stationary near the conformal boundary so as to allow for the 
discussion of gravitational radiation — the existence of this type of spacetimes has been shown 
in [14] . The examples constructed in |14] are somehow special, as they arise as the development 
of Cauchy initial data sets which are exactly Schwarzschildean in the asymptotic end. More 
generally, recent developments in the construction of solutions to the Einstein constraint equations 
by means of gluing methods — see e.g. [13 [TH] HZ] — allow to construct asymptotically simple 
spacetimes from initial data sets which are exactly stationary in the asymptotic region. 

Given this state of affairs, it is natural to ask whether there are more general types of asymp- 
totically simple spacetimes than the ones described in the previous paragraph — see e.g. [37J . This 
question leads to the so-called problem of spatial infinity. If asymptotically simple spacetimes are 
to be constructed using some version of the Cauchy problem in general Relativity, then one has 
to prescribe some initial data set (S, h a b, Xab) on an asymptotically Euclidean hypersurface — for 
simplicity, here it will be assumed that S has the topology of R 3 . The question now is: how does 
one encode in (<S, h a b,Xab) that the development will be asymptotically simple? The examples of 
[13] suggest that this issue has to be related in some way to the behaviour of the initial data set 
in its asymptotic region. 

1.2 Asymptotically Euclidean and regular initial data sets 

As we are discussing properties of spacetimes by means of the conformally rescaled setting given 
by the notion of asymptotic simplicity, it is also natural to work with a conformally rescaled 
one-point compactification of the initial hypersurface S. To this end, we recall the notion of 
asymptotically Euclidean and regular Riemannian manifolds. The pair (S, h a b) will be said to be 
asymptotically Euclidean and regular if there exists a 3-dimensional, orientable, smooth compact 
manifold (<S, h a b), a point i g 5, a diffeomorphism <f> : S \ {i} — > S and a function ft G C 2 (S) n 
C°°(S \ {i}) with the properties 

n(») = o, D„n(i) = o, D a D b n{i) = -2h ab (i), 

Ci>0 onS\{i}, 

hab = Q 2 (t>*h ab - 

Again, in order to simplify the notation, the last condition will be written as h a b — £l 2 h a b so 
that S \ {i} and S are identified. The key property of asymptotically Euclidean and regular 
manifolds is that suitable neighbourhoods of the point i (the point at infinity) are mapped into 
the asymptotic end of S. Thus, one can use local differential geometry to discuss the asymptotic 
properties of the initial data set (iS, h a b). 

1.3 Static and stationary spacetimes 

The notion of asymptotically Euclidean and regular manifolds has been crucial to understand the 
asymptotic properties of static and stationary spacetimes — see e.g. [11] [8] [22] — and to prove 
results concerning their multipolar structure — see e.g. [2 [5] HJ [6] [3] [12] [13l [28]. Stationary (and 
static) spacetimes are best discussed in terms of a quotient manifold X obtained by identifying 
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points on Ai lying on the same orbit of the stationary (static) Killing vector From this 
symmetry reduction procedure one also obtains a metric j a b for the quotient manifold X. As the 
stationary spacetime arises as the development of an asymptotically Euclidean initial data set 
(S, h a b, Xab), the pair {X, j a b) will also be asymptotically Euclidean. Conversely, one can prescribe 
the leading asymptotic behaviour of the initial data set (<S, h a t, x a b) from assumptions on the 
asymptotic behaviour of (X^ajf). For example, one can assume that (X^ai) is asymptotically 
Euclidean and regular and work on a point-compactified manifold X and a conformally rescaled 
quotient metric 7 a f, — this is the assumption made, for example, in 

One of the key results of the theory of stationary spacetimes is that there exists coordinates 
in a suitably small neighbourhood of i for which "f a b, fi and the (rescaled) stationary potentials 
are analytic. This analyticity in a neighbourhood of the point at infinity of the compactificd 
quotient manifold is the key to establish that static and stationary spacetimes are asymptotically 
simple. For static spacetimes, the analyticity on X is inherited by the point compactification S 
of time symmetric slices, and the conformal factor fl of the conformally rescaled quotient metric 
7ab is used as conformal factor for the whole spacetime, so that the spacetimes are asymptotically 
simple. The situation for stationary spacetimes is more delicate: in this case the (analytic) 
conformal quotient metric j a b is no longer conformally related to the conformal metric h a f, of the 
i-constant slices. Furthermore, h a t, is no longer analytic, but only of class C 2 ' a . Notwithstanding, 
it is still possible to construct a smooth conformal extension of the stationary spacetime. This 
result shows that although asymptotic simplicity is a property which can be naturally expected 
from stationary spacetimes, the fact that it holds is a consequence of the structural properties of 
the stationary equations at spatial infinity. 

1.4 The cylinder at spatial infinity 

In order to answer the question of whether there exist asymptotically simple spacetimes arising 
from initial data sets which are neither static nor stationary in a neighbourhood of infinity, one 
needs a framework that allows to resolve and disentangle the delicate structure of spacetime in 
this region. Furthermore, as the strategy to construct asymptotically simple spacetimes is to make 
use of the Cauchy problem in General Relativity, one would like to be able to formulate an initial 
value problem with data prescribed on the compact manifold S for various conformal fields which 
would directly render the conformally rescaled manifold {M.^g llu ). Appropriate tools for this 
construction are the conformal Einstein field equations — see e.g. |19[ 1201 121] — and extensions 
thereof — see [UJ 123 123 HZ]- However, the representation of spatial infinity as suggested by the 
point-compactification of the initial hypersurface S presents us with technical difficulties. The 
underlying reason is that for initial data sets with non-vanishing ADM mass, spatial infinity is 
a singular point of the conformal geometry — see e.g. [22] . At the level of the conformal field 
equations and the various fields they govern, this singular behaviour of the conformal structure 
translates into the divergence of the so-called rescaled Weyl tensor at i. 

In order to overcome the difficulties at spatial infinity that have been described in the previous 
paragraph, a new conformal representation of the region of spacetime near null and spatial infinity 
was introduced in |24) . This representation, based on general properties of conformal structures, 
together with the extended conformal field equations allows to introduce a regular finite initial 
value problem at spatial infinity for which both the equations and their initial data are regular at 
the conformal boundary. 

Whereas the standard (Penrose) compactification of spacetimes considers spatial infinity as 
a point, the approach used in [241 I27j represents spatial infinity as an extended set with the 
topology of [—1,1] x § 2 . This cylinder at spatial infinity is obtained as follows: one starts 
with the standard point-compactification S of an asymptotically Euclidean initial data set <S. 
As in previous paragraphs, S contains a special point i representing the infinity of S. In a 
second stage, the point i is blown up to a 2-sphere. This blowing up is achieved by lifting a 
neighbourhood B of i to the bundle of orthonormal frames with group 0(3) — or equivalently 
to the bundle of space-spinors with group 5[/(2,C). In a final step, one uses a congruence 
of timelike conformal geodesies to obtain a conformal analogue of Gaussian coordinates in a 
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spacetime neighbourhood of B. Timelike conformal geodesies are conformal invariants which 
preserve their quality of being timelike under conformal transformations. Conformal Gaussian 
systems based on these curves provide a canonical conformal class of conformal factors for the 
development of the initial data. Remarkably, these conformal factors can be written entirely 
in terms of initial data quantities. Hence, the location of the conformal boundary is known a 
priori. The conformal boundary rendered by this class of canonical conformal factors contains 
a null infinity with the same structure as in the case of the standard Penrose compactification. 
Spatial infinity, however, extends in the time dimension — so that one can speak of the cylinder 
at spatial infinity proper. Of crucial relevance are the critical sets {±1} x S 2 — the collection of 
points where null and spatial infinity intersect. Null infinity and spatial infinity do not intersect 
tangentially at the critical points. As a consequence, some of the propagation equations implied 
by the conformal field equations degenerate at the critical points. The analysis in [24] — see 
also [4TJ [40[ [42j [43] — has shown that, as a result, the solutions to the conformal field equations 
develop certain types of logarithmic singularities at the critical sets. These singularities form 
an intrinsic part of the conformal structure and cannot be regarded as a consequence of a bad 
gauge choice. The hyperbolic nature of the conformal propagation equations suggests that these 
singularities will propagate along null infinity, and thus, they will have an effect on the regularity 
of the conformal boundary. 

The construction of the cylinder at spatial infinity bears some relation to other approaches in 
the analysis of the structure of spatial infinity. For example, the blow up of the point at infinity 
is closely related to Geroch's idea of directional dependent tensors — see [2U [32 • This idea was 
latter retaken in the discussions given in [2 , 38] . The cylinder at infinity is closely related to the 
hyperboloid of spatial infinity also introduced in [2] , and latter retaken by [9] in a first attempt 
to combine geometric and partial differential equations points of view to study the structure of 
spatial infinity. 

1.5 Static spacetimes and the cylinder at spatial infinity 

In view that static and stationary spacetimes provide prime examples of asymptotically simple 
spacetimes, one also expects the associated construction of the cylinder at spatial infinity to be 
as smooth as it can be. This smoothness can be regarded as a consistency of the setting. If static 
or stationary spacetimes were to exhibit some type of non-smooth behaviour at the cylinder at 
spatial infinity, these by necessity have to be associated to a bad gauge choice. In [27J a proof of 
the smoothness of the cylinder at spatial infinity for static spacetimes was given. Surprisingly, this 
proof is much more complicated than what one would expect given that: firstly, the conformal 
fields are analytic in a neighbourhood of spatial infinity; and secondly, that the spacetimes are 
time independent. The difficulties in the analysis can be explained, in part, by the fact that the 
conformal geodesies used in the construction of the cylinder at spatial infinity are not aligned 
with the orbits of the timelike Killing vector. Nevertheless, the fact that one is considering time 
symmetric spacetimes simplifies the analysis as the quotient manifold can be identified with the 
slices of constant t so that the analyticity of the point i is inherited by spatial infinity i and by 
all spacetime quantities. 

It should be pointed out that the relevance of the analysis of static spacetimes carried out 
in [57] goes beyond its role as a consistency check of the framework. The analysis in [HI [35] 
suggests that static spacetimes have a special position among the class of spacetimes with a 
smooth compactification at spatial infinity. More precisely, it is conjectured that: 

Conjecture 1. // an analytic time symmetric initial data set for the Einstein vacuum equation 
yields a development with a smooth null infinity, then the initial data set is exactly static in the 
neighbourhood of spatial infinity. 

The rigidity results of [HJ H5J constitute further evidence in support of the conjecture. 
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1.6 Our main result 



In the present article we extend the analysis of the cylinder at spatial infinity carried out in |27j 
to the case of of stationary spacetimes. More precisely, we show that: 

Theorem. Given an initial data set for the vacuum Einstein field equations which is stationary 
in a neighbourhood of infinity, the solutions to the regular finite initial value problem for the 
conformal field equations at spatial infinity is smooth in a neighbourhood of the cylinder at spatial 
infinity, and in particular through the critical sets. 

The strategy to prove this result is as follows. Starting with a generic asymptotically fiat 
stationary spacetime (.M, (jr^,,), one makes use of the analysis of [18] to construct a conformal 
completion, (£l, c^ t „), of the stationary spacetime in a neighbourhood of spatial infinity. This 
completion is adapted to the stationarity of the spacetime, and has a smooth null infinity. This 
representation is, however, not suitable for our purposes as spatial infinity is represented as a 
point. The conformal metric g^ is given in terms of some asymptotically Cartesian coordinates. 
In a second stage one performs a change of coordinates to a polar system in which asymptotic 
expansions can be analysed in a more convenient way. A subsequent change of the time coordinate 
and an associated conformal rescaling render a conformal representation, (Q, g^ u ), in which spatial 
infinity already appears as an extended set with the topology of a cylinder. This representation 
is, in a strict sense, not a conformal completion as spatial infinity is not at a finite distance with 
respect to the conformal metric g^ — the metric becomes singular there. In order to deal with 
this behaviour one introduces a suitable frame v a . It is then show that the components of the 
metric, g a f3, with respect to this frame are regular at infinity. Furthermore, the components of 
key derived objects (the Schouten and Weyl tensors) are also shown to be regular. This is the 
most calculational involved part of our argument. Once the cylinder at spatial infinity has been 
obtained, one shows that the cylinder itself, and a neighbourhood of it can be ruled by means 
of a congruence of conformal geodesies. This cannot be shown explicitly, and thus, one has to 
resort to a perturbative argument. This construction leads to the canonical conformal factor 0, 
related to ft through a further conformal factor IT. The congruence of conformal congruences 
also gives rise to a Weyl connection V. An abstract integration of the conformal evolution 
equations along the congruence of conformal geodesies shows that solutions to the conformal 
evolution equations extend smoothly through the cylinder at spatial infinity and also through a 
suitable neighbourhood of null infinity. In a final step, it is then shown that the construction is 
independent of the conformal gauge used to write the stationary initial data. 

Contrasted with the result for static spacetimes given in |27j , the main difficulties in proving 
our main result are: the presence of a non-vanishing second fundamental form in the slices of 
constant t has as a consequence a Weyl tensor with non- vanishing magnetic part; and crucially, 
the quotient manifold cannot be directly identified with the structures of the constant t slices. In 
particular, as already discussed, the analytic structure of the quotient manifold in a neighbour- 
hood of infinity is not inherited by the slices of constant coordinate t. Instead, one obtains fields 
which are of the form f + pg with /, g analytic and p a suitable radial coordinate — recall that the 
radial coordinate is not smooth in a neighbourhood of i with respect to Cartesian coordinates. 
It is of interest to notice that our analysis requires the explicit computation up to quadrupolar 
order of the expansions of the relevant conformal fields. 

Our argument assumes, a priori, the existence of the stationary spacetime, and makes state- 
ments about the smoothness of the spacetime in a certain gauge from the known smoothness 
in another gauge. A proof that makes only use of the conformal evolution and properties of 
stationary initial data sets would be much more complicated and would require a much deeper 
understanding of the properties of the conformal field equations and associated conformal struc- 
tures than the one that is currently available. We expect our analysis to shed some light in this 
direction. 

As in the case of static spacetimes, our main result, on the one hand, ensures that the con- 
struction of the cylinder at spatial infinity for spacetimes without time reflexion symmetry does 
not have spurious gauge singularities, and on the other hand, it is expected to play a key role in 
a proof of a suitable generalisation of Conjecture [TJ 
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Overview of the article 



In Section [5] we present a concise discussion of the Conformal field equations and conformal 
geodesies. The presentation in this section is aimed at providing a context for the analysis of 
the subsequent sections of this article. Section [3] briefly summarises the so-called regular initial 
value problem at spatial infinity. This discussion includes, in particular, the construction of the 
so-called cylinder at spatial infinity. Section[4]discusses results about stationary spacetimes which 
are relevant for our analysis. Particular attention is paid to their asymptotic expansions in both 
the quotient manifold and in a Cauchy slice. A "standard" conformal completion of stationary 
spacetimes is discussed. Section [5] discusses an alternative conformal completion for stationary 
spacetimes. This particular completion ultimately leads to the cylinder at spatial infinity. A 
discussion of the asymptotic expansions for the relevant field quantities in this conformal com- 
pletion are provided. In particular, it is shown that the components of the Schouten and Weyl 
tensors with respect to a particular frame are regular at infinity. Section |B] provides a discussion 
of the construction of conformal Gaussian systems in the neighbourhood of the cylinder at spatial 
infinity of stationary spacetimes. As a result of this analysis, it is shown that in a certain gauge 
the setting of the initial value problem at spatial infinity is as regular as it is to be expected. In 
a second stage it is shown that the construction is independent of the particulars of the choice of 
conformal gauge. This discussion completes the proof of our main result. Section [7] provides some 
concluding remarks concerning our analysis. Some lengthy expansions are presented separately 
from the main text in Appendix [A"l 

2 Conformal field equations and conformal geodesies 

The regular finite initial value problem at spatial infinity, presented in Section [3J was introduced 
by Friedrich in J4| and is based on a conformal representation of the Einstein field equations, 
known as the extended conformal field equations. In this section we elaborate further on the ideas 
discussed in Subsection [O] of the introduction, and we present a concise discussion of this system 
and of its associated structures. The presentation is geared towards the purposes of the present 
article. 

2.1 Weyl connections 

Let M denote a 4-dimensional manifold endowed with a Lorentzian metric g^ u . A conformal 
rescaling of the metric is given by 

where G is a positive function on M.. The conformal class [g] is the collection of all metrics 
conformally related to g^ 

[.<?] = {9^ I 9^ = @ 2 .<?m^ © > o}. 

Let V denote the covariant derivative operator of a torsion- free connection on A4. This connection 
is called a conformal connection or Weyl connection for [g] if for g^ £ [g] one has that 

V^„a = -2 ffj, g v \ (1) 

for some smooth 1-form f^. The connection V preserves the conformal structure of [g\. Further- 
more, it does not depend on the class representative. For example, if g^ v = Q 2 g^ v , then 

with 

/„ = /„- e-%e. 

If V is the Levi-Civita connection of g^, then from equation (Q]) we have that the connections V 
and V define the difference tensor V — V = S(f), given by 

$(/)„ 9 u = 5/ fu + s v " U - g^g pX h, (2) 

where S^ p is the Kronecker delta. This shows that V can be specified using V and / M . 
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2.2 The extended conformal field equations 

We now specialise to the case where g^ u is a solution of Einstein vacuum held equations on A4. , 
Ric[g AI „] = 0. The Weyl connection V and the Levi-Civita connections V and V are related by 

v-v = s(/~), v- v = s(e- 1 de), v-v = $(/). 

The Riemann tensor of the Weyl connection V can be decomposed as 

R M v\ P = 2(5 M [\L p ] v - 5 M v L[x P \ - 9u[\L p ] /1 ) + C M v \ pi 

where L M „ and C 7 ' u \ p denote, respectively, the Schouten and Weyl tensors. The Schouten tensor 
is given by 

where 

Rjili — R^ ppv : R — 9^ Rp,v ■ 

In the sequel it will be convenient to consider the 1-form 

d„ = o.fp = Qf„ + v M e, 

and the rescaled conformal Weyl tensor 

W v x p = 9- 1 C' i v x p . 

In order to deal with the possible direction dependence of the various fields near space-like 
infinity, it is convenient to use a frame formalism and a suitably chosen orthonormal frame field. 
For this, let us take a frame field {e a } a =o,i,2,3 satisfying 

g a p = g(e a ,ep) = rj a p = diag(+l, -1, -1, -1). 

Let V a and V a denote, respectively, the covariant derivatives in the direction of e a with respect 
to the connection V and V. The respective connection coefficients are defined by V a ep = Ta 1 p e 7 
and V a ep = r a 7 ^ e 7 . The frame coefficients with respect to an as yet unspecified coordinate 
system x^ are given by e^ a — (dx^, e a ). Using equation ([2} we have 

f a 1 p = r Q ," l ' ( g + 8a 1 fp + Sp 1 f a - g a /3g' yS fs, 

where f a = a fp,- Then, if all tensor fields are expressed in terms of the frame field and the 
corresponding connection coefficients, the extended conformal field equations are equations for 
the unknowns 

u = (e" a , f a 7 /3, L a p, W a p yS ), 

and are given by 

[e a , ep] = {Ta^p - Tp 1 ' a )e-y, (3a) 

e a (f p S -y) — epipa- 5 -y) — {Fci^p — f p <L a )F e S ~ / + T a S e Tp e y — tp S e t a <L - v 

= '2{g S [cxLpj- f - g S yL[a0] - g~,[ a Lp] s } + QW s ^ a p, (3b) 
V a i/3 7 - VpLa-, — d s W s lot p, (3c) 
V s W\ a p = 0. (3d) 

In the last equation, called the Bianchi equation, one has to consider the relation between To 1 p 
and To 1 p, whence f a = \T a ^ p- Notice that no differential equations are given for the fields 
and d a . This is due to the conformal gauge freedom introduced into Einstein's field equations by 
considering general Weyl connections and conformal metrics. 
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2.3 The conformal Gauss gauge 



The fields and d a can be specified by means of a choice of gauge. To this end we consider 
conformal geodesies. A conformal geodesic for (M,g^) is a curve x(t) in M. and a 1-form /(r) 
along the curve, which solve the following system of ordinary differential equations: 

(V ± ±r + S(f) x %x x x p = 0, (4a) 
(V*/)„ - \U S(/V, x x = L\ u x x . (4b) 

Conformal geodesies are invariants of the conformal structure in the following sense: if x(t), /(t) 
solve equations (|4aj) - (|4bp and b is a smooth 1-form field on Ai, then the pair x(t), (/ — &)U( T ) 
solves the same equations with V replaced by V = V + 5(6) and L by L. This means that x(t), 
and in particular the parameter r, do not depend on the Weyl connection in the conformal class 
which is used to write the conformal geodesic equations. 

Conformal geodesies, and in particular congruences of conformal geodesies, can be used to 
construct a special gauge for the conformal equations. For this, let S be a space-like hypersurface 
in the given vacuum solution (Ad,^^). We choose on S a positive 'conformal factor' 0„, a frame 
field e a *, and a 1-form /*, such that g(e a *,ep*) = Oj 2 ?7 a /3 and eo* is orthogonal to S. Then 
there exists through each point x* £ S a unique conformal geodesic x(t), /(t) with r = on 
S which satisfies there the initial conditions x = eo*, / = /*• These curves define a smooth 
caustic free congruence in a neighbourhood U of S if all data are smooth. Furthermore, / defines 
a smooth 1-form on U which supplies a Weyl connection V on hi given by V = V + S(f). A 
smooth frame field e a and a conformal factor G are then obtained on hi by solving 

V±e a = 0, 

v*e = e<i,/>, 

for given initial conditions e a = e a *, = 0* on S. The frame field is orthonormal for the 
metric — Q 2 g^ u - Dragging along the congruence the local coordinates x a , a = 1,2,3 on S 
and setting x° — r we obtain a coordinate system. In this gauge one has in Li that 

x = e = d T , to^a = 0, L 0a = 0. (5) 

This choice of coordinates, frame field, and conformal gauge will be referred to as conformal 
Gauss system. Remarkably, in this gauge it is possible to obtain explicit expressions for and d 
in terms of the initial data, given by 

9(r) = 8* (l + (/., e *>r + \gHf*J*)r 2 ) , (6a) 

d = 6, d a = 6*(/„ e a *>, a = 1,2, 3, (6b) 

where the quantities with the subscript '*' are constant along the conformal geodesies and <j" 
denotes the contravariant version of g^ LU . These expressions, together with equations (|3a[) - (l3d[) 
provide a complete system of equations for u, called the general conformal field equations. Setting 
a = in (|31i)) - ((3"d"|) and observing the gauge conditions (J5J) one obtains the following evolution 
equations: 

d T e^ a = -f a W (7a) 

Orfa 1 ' (3 = -t S 1 /3f a S + g~* oL a f3 + ' pLo a - g^oP + ® W 1 p 0oil (7b) 

d T L a/3 = d^W 1 0,0(3- (7c) 

If one extracts from the Bianchi equation (|3dp a symmetric hyperbolic system, one gets symmetric 
hyperbolic reduced equations for those components of u which are not determined explicitly by 
the gauge conditions. The resulting system is called the reduced conformal field equations. It can 
be shown that for such a choice of reduced equations, any solution which satisfies (|3a[) - (|3d[) on 
a suitable spacelike hypersurface does indeed satisfy the complete set of field equations in the 
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part of the domain of dependence of the initial data set where is positive. This equations are 
equivalent to Einstein field equations in the sense that if one has a solution of the conformal 
system, then one has a solution of Einstein field equations in the region where the conformal 
factor is positive, and vice versa. The substantial advantage of using the conformal system is 
that one can deal with regions where the conformal factor vanish. Furthermore, in the Gauss 
gauge the location of this region can be prescribed a priori, giving us full control on the conformal 
boundary if the evolution extends far enough. 

3 The regular finite initial value problem at space-like in- 
finity 

As mentioned in the introduction, the construction of the regular finite initial value problem at 
space-like infinity consists of two main steps. For this, one considers a hypersurface S, which is 
a one-point conformal compactification of the Cauchy hypersurface S and therefore contains a 
geometrically distinguished point i — cfr. the discussion in Subsection 11.21 In the first step of 
the construction, the point i is blown up into a sphere TP. In the second step, a congruence of 
conformal geodesies is used to describe the evolution of the conformal fields in a neighbourhood 
of 1° . The construction is described in detail in [27J , and it is implemented through the bundle 
of spin frames over S near i. In what follows, we present a summary of the aspects of a similar 
construction based on the bundle of orthonormal frames — see also [25] . 

3.1 The blow up of spatial infinity 

Start by considering a space-like Cauchy hypersurface S with intrinsic metric h a b containing the 
distinguished point i. Next, choose a fixed oriented /i-orthonormal frame e a , a = 1,2,3, at i. 
Any other such frame at i is obtained by a rotation of e a . That is, all other /i-orthonormal 
frames at i are of the form e a (s) = s b a eb with s = (s b a ) G 50(3). In particular, ea(s) covers 
all possible directions at i as one lets s exhaust 5*0(3). For a given value of s, one distinguishes 
ea(s) as the radial vector at i. Keeping s fixed, we construct the /i-geodesic starting at i that 
has tangent vector es(s) and denote by p the affine parameter on the geodesic that vanishes at i. 
The frame e a (s) is then parallelly transported along the geodesic. For a particular value of the 
affine parameter p, the frame thus obtained will be denoted by e a (p,s). We will consider only 
the region \p\ < a, where a is chosen such that the metric ball B centered at i with radius a is 
a convex normal neighbourhood for the 3-metric h. The map from the set (—a, a) x 50(3) into 
the bundle SO{S) of oriented orthonormal frames over S, given by (p, s) — > e a (p,s), defines a 
smooth embedding of a 4-dimensional manifold into SO(S). 

In what follows, only non- negative values of p will be considered. Denote by B the image of the 
set [0, a) x 50(3). The boundary of B will be denoted by TP. One has that 2° = {p = 0} ~ 50(3). 
Finally, let tt denote the restriction to B of the projection of SO(S) onto S. In the sequel it will be 
convenient to consider the subgroup 5*0(2) of 50(3), given by 50(2) = {s' G 50(3) | s' b 3 e& = e 3 } 
— this is the subgroup of 50(3) whose action leaves es invariant. Accordingly, if s G 50(3) and 
s' G 50(2) then e Q (s) and e a (s s') are parallelly transported along the same geodesic. Hence, 
when we consider the projection n we have that n(e a (p, s)) = n(e a (p, s s')), and therefore the 
map 7r has a factorisation 

B^B' = B/SO{2) ^ B. 

On the one hand, the projection tt" maps the set 7r'(I°) — 5 2 onto i and on the other it implies 
a diffcomorphism of B'\tt'(T°) onto the punctured ball B\{i}. This diffcomorphism can be used 
to identify these sets. However, instead of this, it is convenient to pull back the initial data on 
B to B via tt. Then B becomes the initial manifold, p and s are used as coordinates on it, and 
its boundary T° is a blow up of i. This manifold has an extra dimension when compared to the 
initial hypersurface B. This extra dimension is given by the action of 50(2). Since all fields have 
a well defined transformation behaviour (spin weight) under such action, on the part of B where 
p > vector fields X, c a (p,s), a — 1,2,3, can be prescribed such that X is generated by the 
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action of SO (2) and the vector fields c a (p,s) satisfy T(ir)c a (p, s) = e a (p,s) — see [37]. These 
vector fields allow the introduction of a frame formalism. 

3.2 Implementing the conformal Gauss gauge 

In the second step of the construction, the development of the data is considered. For this, one 
uses a conformal Gauss system — see Section [5] In what follows, it is assumed that the conformal 
compactification of the hypersurface S has been achieved by means of a conformal factor Q. The 
initial data for the conformal factor for the Gauss system is set by requiring that 

e* = k^vl. (8) 

where the function k satisfies 

k = P k', k' e k'>0, Xk' = 0, k'\ i0 = 1. 

The change of the conformal factor induced by the function k implies a map S : e a — > ne a which 
maps the set B\I° bijectively onto a smooth submanifold B* of the bundle of frame fields over 
B. The diffeomorphism S is used to carry the coordinates p, s and the vector fields X, c a to B*. 
The projection of B* onto B will be denoted again by it. 

The construction of the conformal Gaussian system requires initial data for the 1-form /. For 
this one takes 

(f,dr)=0, 7T*f = K- 1 dK. (9) 

The reduced field equations (|Ta[) - (|7c[) and the symmetric hyperbolic system implied by (I3dp 
can be interpreted as equations in the development of B*. This development will be denoted by 
TV, and is a 5-dimensional manifold smoothly embedded in the bundle of fame fields over M. 
The manifold TV" is a 50(2) bundle over the spacetime. Its projection sending TV" onto M. will 
again be denoted by 7r. The coordinates p, s and the vector fields X, c a are pushed forward with 
the flow of the conformal geodesies ruling TV, in such a way that X generates the Kernel of ir. 
The parameter X° = T defines a further independent coordinate with x° = r = on B*. The 
tangent vector field of this congruence is denoted by d T . To interpret the reduced field equations 
as equations on TV it is assumed that frame fields e a are vector fields on TV which are defined at 
a frame field (d T ,c a ) by the requirements that: (i) they project onto the frame field defined by 
(<3 T ,c a ) on Ai; (ii) they do not pick up components in the X direction. The unknowns in the 
reduced field equations are then interpreted as vector- valued functions on TV. 

3.3 The conformal boundary 

From equations (j5a]) and ([§]) it follows that 

e = e*(i-r 2 ^, on tV, (io) 

where the function u> is given by 

ui = — =, on B . 

y/\DJWHT\ 

As before, subscripts '*' imply that the relevant functions are constant along the conformal 
geodesies. Explicit expressions for d a can be obtained from (|6bj) using that / and / are related 

by/ = /-e- 1 de. 

An important property of the construction described in the previous lines is that if the initial 
data set for the reduced conformal field equations has a smooth limit as p — > 0, then it can be 
smoothly extended into the coordinate range p < 0. Similarly, and d a take smooth limits 
as p — >• and can then be extended smoothly into a range where p < 0. It follows that the 
initial value problem for the reduced field equations can be extended smoothly into a range where 
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p < in such a way that the reduced equations are still a symmetric hyperbolic system. If the 
development of the initial value problem just formulated extends far enough, then the following 
regions of the development can be distinguished: 

Af= ||r| < -, < p < a, s e 50(3)} , 

Af= {|t| < -, < p < a, s e 5(9(3)} , 

where w/k is a function of p and s. One also has the 4-dimensional submanifolds 

= {r = +-, < p < a, s G 50(3)} , 

Jf" = |t = --, < p < a, s e 50(3)} , 

2= {t < 1, p = 0, s e 50(3)}, 

and the 3-dimensional submanifolds 

1+ = {t = +1, p = Q, se50(3)}, 
I - = {r = -1, p=0,se 50(3)} , 
1° = {t = 0, p = 0, s e 50(3)} , 

where it has been observed that ui/k — >lasp— > 0. It can be verified that 

6 > on TV", 

6 = 0, d6 / on/^uZ+Ul, 
6 = 0, d6 = onrul + . 

The initial hypersurface B* is given by 

B* = {t = 0, < p < a, s G 50(3)} . 

Its closure in A/" is given by 

B = {t = 0, < p < a, s e 50(3)} = B* U Z°. 
Factoring out the group 50(2) projects N onto the set M., representing the "physical spacetime" . 

Observation. The data on B* have a unique smooth extension to B. Furthermore, the solution 
on N depends only on this data as the set I is a total characteristic of the system of equations, 
and therefore the solution there depends only on the data on 1°. The set X, referred to as the 
cylinder at spacelike infinity, represents a boundary of the spacetime TV, and may be understood 
as a blow up of spacelike infinity i . Of particular importance are the sets X + and X~ , the critical 
sets, as the system of evolution equations degenerates there. This degeneracy makes very difficult 
to make any statement about smoothness of the solution to the initial value problem stated above, 
even it has been shown that for general initial data this sets are not regular — see e.g. [2~il |2"T] . 

In view of the discussion of the previous paragraphs, the objective of the present article is as 
follows: to show that for stationary asymptotically flat initial data the solutions to the regular 
finite initial value problem are smooth in a neighbourhood of Z, and in particular are smooth 
through 1 + and I~ . 

4 Stationary asymptotically flat spacetimes 

In what follows, let (M, g^, £ M ) denote a a stationary spacetime. That is, Ai is a four-dimensional 
manifold, g^ v is a Lorentzian metric on M. with signature (+,—,—,—) satisfying the Einstein 
vacuum field equations, and £ M is a time-like Killing vector field with complete orbits. As discussed 
in Section 1 1 .31 instead of working with the 4-dimensional manifold M. it is more convenient to 
consider the quotient manifold, X, of M with respect to the trajectories of the Killing vector £^ 
—see Ha 123 ED]- 
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4.1 The stationary metric in terms of quantities on the quotient man- 
ifold 

Locally, the metric g^ v can be written in terms of quantities denned on the quotient manifold X: 
a scalar V, a 1-form f3 a , and a Riemannian metric ^f ab . More precisely, one has that 

g^dx^dx" = V(dt + p a dx a )(dt + p~ b dx b ) - V- 1 % b dx a dx b , (11) 

where V, P a and ^f ab depend only on the spatial coordinates x a . 

In order to obtain the field equations implied on X by Einstein vacuum field equations on M. 
it is convenient to consider the quantity Cj a , defined on X by 

Cj a = -V 2 e abc D b p~ c , (12) 

where D is the covariant derivative associated with "f ab . Then the Einstein vacuum field equations 
on M. imply 

D[ a u>b] = 0. 

If one further assumes X to be simply connected (our analysis will concentrate in a neighbourhood 
of infinity), then there exists a scalar field uj such that 

D a uj = Q a . 

In the sequel we will consider 'gauge' transformations of the form 

Pa -> Pa + a f, 

where / is a scalar field on X. Clearly, oj a does not change under these transformations — cfr. 
equation (fT2]) . Moreover, the metric remains unchanged if one sets t — > t — f. 

4.2 The Hansen potentials 

In order to write down the stationary field equations it is convenient to introduce the so-called 
Hansen potentials: 

7 V 2 +uj 2 -1 - w ~ V 2 +uj 2 + 1 



} M = 7T7 , PS = 77F7> 9K 



4V ' r 2V" W 

They are not independent as 

4>m + 4>s - 4>k = -i 

The vacuum field equations on M. then imply on X 

A0 Q = 2R[^ a , a = M, S, K, (13a) 

R ab [j] = 2{D a j> M Db4>M + Da4>sD b 4> s ~ D a ^ K b h ~ VK ). (13b) 

The latter will be regarded as field equations for 7 a f,, (f>M and <ps on X. They are equivalent to 
Einstein vacuum field equations on M in the sense that A4 can be reconstructed as a stationary 
spacetime if 7 a f,, 4>m and 4>s are given. 



4.3 The 3+1 form of the stationary metric 

As mentioned in Section ri.3l of the introduction, although the field equations take the simple form 
(I13a[) - (|13b>[) in X, our main interest is to consider the Cauchy problem with data arising from 
stationary spacetimes. For this one needs to consider the 3 + 1 decomposition of the spacetime 
metric with respect to a particular spacelike hypersurface of Ai. If we choose S as defined 
by t — constant, then g^ v has a 3 + 1 decomposition with respect to S given by 

g^dx^dx" = N 2 dt 2 - h ab (N a dt + dx a )(N b dt + dx b ), (14) 
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where N, N a , h a b denote, respectively, the lapse function, the shift vector and the intrinsic metric 
of the hypersurface S. Comparing the line elements (TTTT) and (fl4|) one finds that 

* 2 =r^sp *"-i J 9|p «-- ^-"AA- < 15 > 

We will adopt the convention of moving indices of objects in X with the quotient metric j a b- 
The relations (1151) allow us to go back and forth between quantities defined on X and quantities 
defined on S. 



4.4 Asymptotic flatness and its consequences 

Following the usual assumptions about asymptotic flatness for stationary spacetimes, it will be 
assumed that (X,"f a b,4>M,4 > s) is asymptotically Euclidean and regular in the sense discussed in 
Section rOl of the introduction — see e.g. jTSJ H3] ■ More precisely, it is assumed that there exists 
a manifold X, such that X = X U {i}, where i is a point. Furthermore, it is assumed that for 
some real constant B 2 > the conformal factor 

ft = ±£T 2 [(l+ 44> 2 M +4^)3-1] (16) 

is C 2 ' a on X and satisfies 

Cl(i) = 0, D a n(i) = 0. 

In addition it will be assumed that 

7afc = £l 2 7afc (17) 

extends to a C 4,a metric on X and satisfies 



D a D b n(i) = 2y ab (i), 

where D is the Levi-Civita covariant derivative of the 3-metric "f a b- 

The conformal rescaling of the metric given by (|17p suggests the following definition of rescaled 
potentials: 

l _ 

(j) a = n ■2 ( p aj a = M,S,K. 

The motivation behind the introduction of conformally rescaled fields is the following theorem by 
Beig & Simon [12] — see also [34] . 

Theorem 1 (Theorem 1 of 12 ). For any asymptotically flat solution (jab, </>m, <As) of the sta- 
tionary equations (|13a[) - (|13b[) there exists a chart defined in some neighbourhood of i in X such 
that (~fab,<ftAi,4 > S,ty are analytic. 



Remark. Given the chart indicated by the previous theorem, one can make a coordinate trans- 
formation to 7-normal coordinates x a centered at i. The fields ("fab, 0a/, 0s, ^) are also analytic 
with respect to the normal coordinates x a — this follows from the Cauchy-Kovalewska theorem 
applied to the equations of the radial geodesies written in the analytic coordinates given by 
Theorem [T] It is important to notice that Theorem [T] does not make any assertion about the 
smoothness of other quantities on X, like V, 4>k, Pa or quantities defined on a hypersurface of 
the spacetime. An analysis of the regularity of these and other related quantities has been carried 
out in PJ]. 

In the sequel, we will require several results from [15] . These will be presented here for 
completeness and quick reference. In the following let the radial coordinate p be defined by 

/ 3 \V2 

In |18] it was found that the non-analyticity of the relevant functions is of a very special type 
and it depends only on the coordinate p. Accordingly, one defines the following function space: 
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Definition (Definition 2.2 of [H]). We define the space E u as the set 

S w = {/ = / 1 +p/ 2 :/ 1 ,/ 2 eC-}, 
where C u denotes the set of analytic functions in a neighbourhood of i. 

Associated to the latter definition one has the following: 
Lemma 1 (Lemma 2.3 of [15]). Let f,g G E u , then 

(i) f + geE". 

(ii) fgeEP. 

(Hi) Iff^O then 1/f G E" . 

Obviously, if / G C u then / G E w . The main result of the analysis in [18] is that most of the 
relevant quantities belong to E" . In particular, one has that: 

Lemma 2 (Lemmas 2.4 and 2.5 of |18|). In the normal coordinates implied by Theorem]]] one 
has that 

V G E". 

Furthermore, there exist a choice of gauge for which the 1-form (5 a has the following form: 

P* = ft + —, (19) 
P 

where j3\, /3 2 are analytic functions of x a given by 

Pl = e abc flx c , (3 2 a =e abc f b x c , (20) 

where /f , f% are analytic and e a bc is the flat volume element. In particular, it follows that 

o. 

4.5 A first conformal compactification of stationary spacetimes 

As discussed in [18] , the fields fi and V can be used to construct a first conformal compactification 
of the physical manifolds M. and S. For this one introduces a conformal factor Ct via 

h = v l ' 2 n. 

Note that, as defined, Ct is not analytic as V is not analytic. The associated rescaled metrics are 
then given by 

9nv = Cl g^v, h a b — Q h a b- 

Hence, one has that 

S^dx^dx" = V 2 n 2 (dt + (3 a dx a ){dt + f3 b dx b ) - j ab dx a dx b . (21) 

Moreover, one also has a 3+1 decomposition with respect to the hypersurfaces S = {t = constant}, 

g^dx^dx" = N 2 dt 2 - h ab {N a dt + dx a )(N b dt + dx b ). (22) 

By comparison with equation fj L4[) one gets that 

N = QN, N a = N a , (23a) 

Kb - lab - V 2 n 2 f3 a (3 b , (23b) 

where h ab is the intrinsic metric on S. Fundamental for our subsequent analysis is the following 
result. 
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Theorem 2 (Theorem 2.6 of 15]). Assume /3„ is given by Lemma\^ Then, in some neighbourhood 
of i, the metric h a b has the form 

kb = h 1 ab + p 3 h 2 ab , (24) 

where h\ b , h ab £ C u . 

Remark. The later result implies that the conformal 3- metric h ab is in C 2 ' a . Therefore, the 
conformal factor Cl can be used to define a conformal compactification S of the Cauchy slice S 
plus the point at infinity i, in the same way as made for X. This implies that the pair (5, h ab ) 
admits a C 2,a compactification — that is, the pair is asymptotically Euclidean and regular in the 
sense discussed in Section 11.21 of the Introduction. The decomposition of f3 a given by equation 
(flU|) in Lemma [U and thus, the decomposition of h given by equation is preserved under the 
transformation 

If one imposes the condition x a /3 a — 0, then one fixes d a f ■ Accordingly, j3 a , as given by equation 
(flUf is the unique possible choice of the 1-form f3 a that satisfies this condition. 

Let x ab , Xab denote, respectively, the extrinsic curvatures of S with respect to the metrics g^ v 
and One has that \ ab — Cl~ 1 x a b- Furthermore, we define 

Ipab = Q^Xab = &~ 2 Xab- 

The behaviour of ip ab near i is given by the following result. 

Theorem 3 (Theorem 2.7 of [18). Assume f3 a as given by Lemma\^ Then in some neighbourhood 
of i, the tensor tp a b has the form 

i>ab = p~ 5 fx {a f3l } + p~ 3 i>l b , 

where ^ ab G , f e E u and /3 2 is given by $2U§. Furthermore, p 8r ipab^ ab e E u . 



4.6 Detailed expansions at infinity 

Key for our present analysis is that if a quantity belongs to the space E^ , then although it is not 
analytic, it nevertheless has in a neighbourhood of infinity an analytic expansion in terms of the 
radial coordinate p and the angular coordinates. In the sequel it will be necessary not only to 
know that relevant fields belong to E u , but also to know the first orders of the expansion in a 
neighbourhood of i. In 37] the first orders of the asymptotic expansions of the unrescaled fields 
4>Mi <Ps, 4>k and jab has been explicitly given in terms of constant tensors M, S, M a , S a , M ab , 
S ab , etc. These expansions read 

M M a x a M ab x a x b M(M 2 + S 2 ) 4 
6 M = — + -Z3- + 2f8 + ~ 3 + 0(r ), 



r r 



- S S a i a S ab x a x b S(M 2 + S 2 ) 4 
4>s = - + -T3- + + + 0{r ), 



1 M 2 + S 2 2MM a x a 2SS a x a 



•>k = - 



2 f 2 f 4 f 4 



+ 0(f- 4 ), 



M 2 _ 2 . 2MM {a x b) 2MM c x c S ab 4MM c x c x a x b 

lab = S ab ~ —(SabT - X a X b ) - - + - 

S 2 _ 2 _ _ . 2SSr a X b \ , 2SS c X C S a b 4,SS c X C X a X b , ~.f~-i\ 

+^(S ab r +x a x b ) + £ + ^ ^ + 0(r ), 

where indices in the coordinates and constant tensors are moved with the flat metric 5 a b- 

Remark. Asymptotic flatness implies that the angular momentum monopole S has to vanish. 
Furthermore, by a suitable choice of the origin for the coordinates x a one can make M a = 0. In 
order to simplify our computations we assume this choice of coordinates. 
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Following the discussion in the previous subsections, we want to consider the fields in a 
neighbourhood of i written in terms of normal coordinates centered at i. This is done in several 
steps. First, one changes coordinates from x a to y a , via 

ia = -fl> r 2 = 6 ab y a y b = ± 
Then, using the expression (IT5|) for the conformal factor with B = M, one finds that 

A = r 2 + M 2 r 4 + ^(S a y a ) 2 r 2 + ±M ab y a y b r 2 + 0(r 5 ). 

Furthermore, from the definition of the conformal potentials if follows that 
<Pm = M + ~M 3 r 2 - ^(S a y a ) 2 + 0(r 3 ), 
<fe = S a y a + \s ab y a y b + 0(r 3 ), 

+« = b + \ M2r - *k W)2 - m- r M ^ ayb + 0{r2) - 

Next, we give the expansion for the conformally rescaled metric jab- However, instead of giving 
these expansions with respect to the coordinates y a , we use normal coordinates, x a , centred at i 
— as it was done in the remark after Theorem [T] For this one requires that after the coordinate 
transformation y a x a the metric satisfies 

j ab x b = 5 ab x b . 

A lengthy computation renders 

lab = S ab + ^ ( M Va& - e a e b ) + 2M (M ab + 6 ab M cd e c e d - 2e {a M b)c e c ) 
+2(S a S b + 5 ab (S c e c ) 2 - 2e (a S b) S c e c )) + 0(p 3 ), 

where p is defined by (fT5)l and e a = x a / p. The leading terms of the expansion of the inverse 
metric are given by 

l ah = <5 Qb - TTTTo ( M 4 (S ab - e a e b ) + 2M(M ab + S ab M cd e c e d - 2e ( - a M b ^ c e c ) 
+2(S a S b + S ab (S c e c ) 2 - 2e (a S b) S c e c )^j + 0(p 3 ), 
The transformation between the coordinates y a and x a is given by 

+ ^pP 3 ( " AfV + M(M ab e b - 2e a M bc e b e c ) 

+S a S b e b - 2e a (S b e b ) 2 ) + 0(p 4 ), (25a) 
r = p- ^p 3 (A/ 4 + MM ab e a e b + (5 a e Q ) 2 ) + 0{p A ). (25b) 
Using (|25a[) - (|25b[) , one can express the rescaled potentials in normal coordinates. One finds that 

4>M =M+ ^P 2 (M 4 - (S a e a f) + 0(p 3 ), 
<fe = ~S a e a p + X -p 2 S ab e a e b + 0(p 3 ), 



y a = x a 
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9 U = ~ - — + ^r7^(H^ 4 - MM ab e a e b - (S a e Q ) 2 ) + O^ 1 ), 



For later use it is also convenient to calculate the expansion of other quantities. Namely, 
V = l + 2Mp + 2M 2 p 2 + -i-p 3 (4M 4 + MM ab e a e b - 2(S a e a ) 2 ) + 0(p 4 ), 

olVl 

f3 a = e abc e c (2S b + 0(p)). 

Using the formula (|21[) one finds that the 4-dimensional spacetime metric and its inverse are given 
by 

r s_( v 2 n 2 v 2 n 2 p a 

M " ^ V 2tfp b - lab + V 2 n 2 f3 a f3 b 

GT) = v 2 ^ 2 p 

\ f3 l -7 ab 

where f3 a = j ab /3 b . Moreover, one has the following expansions for the components of the metric 
and its inverse: 

g tt =p 4 + AMp 5 + JL p 6 (l3M 4 + MM ab e a e b + (S a e a ) 2 ) + 0(p 7 ), 
gta = P i e abc e c (2S b + O(p)), 

9ab = -S ab - -£p( Mi ( 6 *b - e a e b ) + 2M(M ab + 5 ab M cd e c e d - 2e {a M b)c e c ) 

+2(S a S b + S ab (S c e c ) 2 - 2e {a S b) S c e c )) + 0(p 3 ), 

1 AM 2 
~^~~p^ + 3M2p 2 
g ta = e abc e c (2S b + 0(p)), 

g ab = -5 ab + -£-^^M i {S ab - e a e b ) + 2M(M ab + 5 ab M cd e c e d - 2e {a M b '> c e c ) 
+2(S a S b + S ab {S c e c ) 2 - 2e (a S ,b) S* c e c )) + 0(p 3 ). 

One also finds that 

fi = V 1 ' 2 ^ = p 2 + Mp 3 + ^M 4 + MM ab e a e b + (S a e a ) 2 ^j + 0(p 5 ). 

5 Conformal extension of stationary vacuum spacetimes 

In this section we introduce a conformal extension of vacuum stationary spacetimes which is 
well adapted for the analysis of the structure of spatial infinity. To this end, we start from 
the conformal extension of the stationary metric given by the line element (|21[) . In order to 
analyse the regularity of the relevant fields we will make use of the first orders expansions in a 
neighbourhood of i, which we collect as we consider the corresponding quantities. The present 
analysis is based on a similar analysis for static spacetimes given in |27j . 

Recall that x a , a = 1, 2, 3, are normal coordinates of the quotient metric 7 b. We have defined 
/ 3 \ 1/2 

For constant t, the surfaces of constant p are diffcomorphic to a 2-dimensional sphere. Accordingly, 
arbitrary coordinates ifj A , A = 2, 3 on the 2-sphere § 2 = {\x\ — 1} can be used to parametrize e a . 
We then write 

e Q = e Q (V A ), de Q = e a dip A . 
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The coordinates ip A can be chosen such that e° depends analytically on them. Consistent with 
the previous definitions one has that x a = pe a (ip A ). Therefore, the conformal 3- metric h a b given 
by equation (|23b|) takes the form 

h = dp 2 + p 2 k, 

where k denotes the 2-dimensional metric on the surfaces of constant p given by 

k = k AB d4i A d^ B = h ab {p)de a de b . 

Remark. Notice that as p — > 0, the metric h approaches the standard Euclidean metric in 
normal coordinates and the metric k approaches the standard line element da 2 = fc(0, ip A ) on the 
2-dimensional unit sphere in the coordinates ip A . 

5.1 Coordinates for the analysis of the cylinder at spatial infinity 

The coordinates (t,p,tp A ) are well adapted to the description of spatial infinity as a point. In 
order to resolve the structure of the cylinder at spatial infinity, it is convenient to introduce 
new coordinates (? 7 p 7 ip A ). The coordinate change is inspired by an analysis of the Minkowski 
spacetime — see e.g. [231 I3H] ■ 

We define x° — t, x° — f, x 1 — p, x A — ip A and consider the map $ : x^ — > x fJ '(x fi ) given 

by 

t(x"') = x°(x"') 1 



, (1 _ f)p {Vn)[se^ A )Y 
x a (x"') = (l-f)pe a [^% 

where the squared brackets indicate the arguments of a given function. One explicitly finds that 

dt = -, TTT-^ — ; rdf+f- — \ — AdS+l, (26a) 

{Vn)[(l -?)pe a ] y(Vn)[pe a ] (Vn)[(l-f)pe a ]J H ' v ' 

dx a = -pe a df + (1 - f)e a dp + (1 - f)pde a , (26b) 

where 

fP ( ^ 



/ - l A di> A , l A = ds, de a = {e a )^ A . 

J(i-f)p \(vn)[se *]J^ A 

The differentials dt, dx a are independent for < f < 1 and p between zero and a small enough 
number. Therefore one can consider the x M as smooth coordinates on an open neighbourhood of 
space-like infinity in {t > 0}. For later use we notice the relation 

p= (1 - f)p, 

which will be used to simplify the notation. 

The main purpose of the coordinate transformation is to remove the coefficient V 2 VL 2 from 
the time-time component of the metric g M „ and to introduce a convenient parametrisation set at 
which p = 0, so that it seems to have an extension in the time direction. 

5.2 A frame adapted to spatial infinity 

We define a set of frame fields v a and their associated coframe fields a a by 

v = d f , fi = pdp, v A = d^,A, 

a = df , a 1 = idp, a A = dip A . 
P 
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To change the frame field associated to our original coordinates to the frame v a one needs the 
inner products v fi a = (dx M , v a ). From the expressions (|26al) and (I26bp one sees that: 



(Vn)[pe a ] 

= 4 - 2M P - + ^7jp(5M 4 - MM ab e a e b - (S a e a ) 2 ) + 0(p 2 ), 
p z p 6M Z 

t = P P 

v 1 {Vn)[pe a } {Vn)[pe a } 

= ~ + ^p f p( 5Mi - MM -be a e b - (S a e a ) 2 ) + 0(p 2 ), 

V 1 A = Ia, 

v a o = -pe\ 
v a i=pe a , 
v a A = pe a ^ A . 



The frame and coframe fields distort the length of the radial component of the tensorial fields 
they are contracted with. This distortion will be of importance in the sequel when discussing 
objects that are singular at spatial infinity. 

5.3 A conformal metric containing the cylinder at spatial infinity 

Let g^u denote a metric conformal to g M „ and g^ u defined by 



It follows that 



v 2 n 2 / P 2 v 2 n 2 (3 a / P 2 
v 2 n 2 p b / P 2 (^ ab + v 2 n 2 p a p b )/p 2 



and 



3M< 



9tt =P 2 + 4Mp 3 
9ta = p 2 e abc (2S b + 0(p))e c , 

9ab 



p 4 (!3Af 4 + MM ab e a e b + (S a e a ) 2 ) + 0(p 5 ), 



P 



1 s 1 

2 ah 3M 2 



M 4 {5ab - e a e b ) + 2M(M ab + S ab M cd e c e d - 2e (a M b)c e c ) 
+2(S a S b + 5 ab (S c e c ) 2 - 2e (a S b) S c e c )\ + O(p). 
The inverse metric g^ v is given by 

(sn - 



P 2 (i - v 2 n 2 (3 c /3 c )/v 2 n 2 P 2 /3 a 

P 2 f3 b 



p 2 J ab 1 ' 



and 



4> - — + ^2 ( 11M4 " MM abe a e b - (S a e a ) 2 ) + 0(p°), 



p" p 3M 2 

ta = p 2 e abc {2Sb + (p))e c , 
4 r 

M 4 (S ab - e a e b ) + 2M(M ab + 8 ab M cd e c e d - 2e {a M b ^ c e c ) 



-ab 



_ p 2 S ab , P 



3M 2 



+2(S a S b + S ab (S c e c ) 2 - 2eS a S b) S c e c ) + 0{p 5 ). 
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Notice that the metric g^ v is singular at p = 0. Thus, the points for which p = are at an infinite 
distance with respect to this metric — hence one does not obtain a finite representation of spatial 
infinity. However, this singular behaviour is counteracted by the use of components with respect 
to the frame and coframe basis introduced in the previous subsection. 

It is also noticed that the conformal factor has the following expansion: 

fl = -V l t 2 tt = p + Mp 2 + -L-p 3 (5M 4 + 2MM ab e a e b + 2{S a e a ) 2 ) + 0{p 4 ). 

p OM 

In the sequel, we will also require the components of the metric g with respect to the frame v a : 

g a p = (®*(g);v a ,vp) 

= ((g^ o ^)dx fl dx v ;v a ,Vp) 

= (g tt °®)(dt,v a )(dt,vp) + 2(g ta o <I>)(dt,v a )(dx a ,Vf3) 

+ {gab ° $)(dx a , V a ){dx b 7 Vp) 

= (gtt ° aV/3 + 2(g ta o $)w* a w Q ^ + {g ab o $)v a a v b p. 
These components are explicitly given by 
.9oo = 0, 
.9oi = 



[Vn)[pe a ] 



goA 



5n 



(l-r) 2 (Vft)[pe°] 
(VQ)[pe a 



(1 -tYp 

[Vtt)[pe a ] ( (V{l)[pe a ] 



giA = 



gAB 



(1 - f) 2 (VCl)[pe a ] \ (Vn)[pe a 
[Vn)[pe a ] ( {VVL)[pe a 



2(1 -r) , 



(1 -f) 2 p \ (VVt)[pe a 
{VVt) 2 [pe a 



(1-f) \{Ia + pPa), 



(IaIb + ZpPaIb) + k AB , 



(i - f ) v 

where j3 A = fl a e a ^A. 

The expansions discussed in the previous paragraphs imply that 

gap = g* a p + 0(p 2 ), 

with 

/ 1 - 2Mp? 

(g a p)=[ l-2Mp? -(l-f)(l+f-4Mpf 2 ) 

V o o k AB (o) 

In order to obtain these expressions we have used that k A B(0) — S a be a ^A e b 1 pB, S a i,e a e b = 1 
and S a b e a e b ^a = 0. The inverse metric is given by 

g a0 = g* al3 + o(p 2 ), 



with 



(g* a>3 ) 



/ (l-f)(l + f-4Mpf 2 ) 1 

(1 - 2Mp?) 2 1 - 2Mpf 



V 



1 - 2Mp? 







„AB 



(0) / 



Remark 1. One sees that although the conformal metric g^ is singular at p = 0, its components 
measured with respect to the frame of Subsection 15.21 are regular and indicate the existence of 



20 



an extended set with the topology of a cylinder at spatial infinity — its sections of constant r 
correspond to 2-spheres. 

Remark 2. The stationary metric g possesses the Killing vector £ = dt- As the conformal factor 
n and p do not depend on t, £ is also a Killing vector of g. In the new coordinates it takes the 
form 

£ = «l ((1 _ f)a , + ^) 
= 2ffiM ((1 _ f) „ 0+ „ l) . 

Remark 3. The metric g is a conformal representation of the metric g which allows us to 
construct an extension of M in a neighbourhood of spatial infinity. For this, we replace the 
hypersurface S by the manifold with boundary S, by adding to S the 2-dimensional surface dS. 
The points of this 2-dimensional surface are thought of as ideal end points of radial curves on S 
as p — s> 0. The coordinates p, ip A extend by definition to smooth coordinates on S and on dS 
we have p = 0. The coordinates ip A , although not specified, are supposed to cover the § 2 . The 
construction described in this paragraph provides an alternative implementation of the blow up 
of spatial infinity discussed in Subsection 11.41 of the Introduction. 

Working by analogy with the discussion of Section [3l one defines the following regions in terms 
of the range of the coordinates f , p and ip A : 

M' = {0 < f < 1,0 < p}, 
M' = {0 < f < 1,0 < p}, 
,J?+' = {? = l,p>0}, 
1' = {0 <r < l,p = 0}, 
Z+' = {r = l,p = 0}, 
X ' = dS = {? = 0, p = 0} , 
X' =l'Ul + '. 

The names for these sets have been chosen in accordance with the related sets defined in 
Section |31 although they differ in some aspects. One can readily verify that in terms of the new 
coordinates, frame and coframe fields, the metric g extends smoothly through the sets and 
I' . Accordingly, M' provides a suitable extension of M.. In order to use this extension to show 
that the construction of the cylinder at infinity for stationary spacetimes is as smooth as expected 
we need also information regarding the Schouten tensor and the conformal Weyl tensor, which 
we derive in the following subsections. Finally, notice that in terms of the coordinates (f, p, tp A ), 
null infinity appears to be parallel to the surfaces of constant f , and in particular the initial 
hypersurface S. 

5.4 Expansions of the Schouten tensor 

In this section we discuss expansions of the Schouten tensor of the metric g. This is related to 
the Schouten tensor of the metric g by 

= L^u - i V M V„0 + V A tt V A ^ g„ v - 

In what follows, we will consider the components of L in the frame v a : 

L a p = ($*(L);v a ,Vp) 

= (L tt o $)w* a w t /3 + 2(L ta o $)w t a w a (3 + (L ab o <$>)v a a v b p . 
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As g is a solution to the vacuum Einstein field equation, it follows that L^ v — 0. Furthermore, 
in view that £1 does not depend on t one obtains 

Alternatively, one can write 

Lu = if t a t d a n + JL § ab d a n d b n g tt , 



L ta = 7 rr t b ad b n + -=-g bc d b nd c ng ta , 

L ab = -^(d a d b n -T a c b d c n) + dc n dd n g ab _ 

Using the expansions for the conformal factor, the first and second derivatives of the conformal 
factor and the Christoffel symbols of the conformal metric given in Appendix [5] ones get that 

L tt = ~p 2 +4Mp 3 + 0(p% 



L ta = -2p 2 e abc e b S c + 0{p 3 ) 

1 ,r „ ^ M 
z-^{o ab - 2e a e b j H 

2p P 



L <*b = TT^i^ah - 2e Q e b ) + — (S a b - 3e a e b ) + O(p ). 



One then has all the ingredients to calculate the components L a p. One obtains 
L 00 = O(p 2 ), 

L 01 = ^+M(2-3f)p + 0(p 2 ), 
Loa = 0(p 2 ), 

in = - f 2 ) - 2M(1 + 2f 2 )(l - f)p + 0{p 2 ), 

Lia = 0(p 2 ), 

Lab = \ k AB + M{\ -f)pk AB + 0{p 2 ), 
where it has been used used that 

S ab e a ^ A e h ^ B = k AB {0) + 0(p 2 ), S ab e a e b ^ A = 0. 
Summarising, one has proven the following result: 

Lemma 3. The components, L a p, of the Schouten tensor in the frame v a are regular (i.e. 
non-divergent) at I'. 

5.5 The Weyl tensor 

In this section we verify that conformal Weyl tensor has a smooth limit at p — » 0. For this, we 
start by recalling the decomposition of the Weyl tensor in terms of its electric and magnetic parts 
with respect to the normal to a hypersurface. 

5.5.1 The electric-magnetic decomposition 

Let M. be the spacetime with metric g^ u and S a space- like hypersurface with unit normal h^. 
The induced metric on S by g^ v is given by 

hfii/ 9^ ri^riy. 

For convenience one defines 

P[il> — n^n U: e^ v \ = nPep^ v \. (^7) 
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The n-electric and n-magnetic parts of the conformal Weyl tensor are given, respectively, by 

Cup = C '^vXpfi^ h x , c* p = C^^ph^h , (28) 

where C* u>i denotes the dual of the conformal Weyl tensor: C* uXp = \Cp Uax e' JX Ap-The n-electric 
and n-magnetic parts of the Weyl tensor are symmetric, trace-free and spatial: 

n%„ = 0, h^ v = 0. 
Given two tensors of rank two, f pu and k pu , their Kulkarni-Nomizu product is defined as 

(/ k)^ Xp = ^{fp[\k p ]y - fv[\k p ]n)- 
In terms of the Kulkarni-Nomizu product, the conformal Weyl tensor is given by 
Cp V \ p = ^{p v [\c p \p - Vp[\c p \ v - n[\c* p]s e 5 pv - n^c*^ s e s Xp ) 
= (P0 c'Tp.uXp - (p c)p uXp - 

5.5.2 Conformal rescalings 

If 9p,v is a solution of the Einstein vacuum field equations, then the first and second fundamental 
forms of S induced by g pv satisfy the Gauss and the Codazzi equations. These equations allow 
to write the pull-back of c pil and to S in terms of the initial data quantities. More precisely, 

Cab = -r ab [h] +Xc C Xab ~XcaXb c , c* ab = ~ D cXd(ah) cd • (29) 

If 9pv = & 2 gpv, then it is well known that v \ p — u \ p . The rescaled conformal Weyl tensor 
is given by v \ p — £l~ 1 C^ v \ p , and therefore 

Wp V \p — QCpvXp. 

The n-electric and n-magnetic parts of W are defined in accordance with (j28|) as 

w vp = Wp„ Xp n>*n x , w* vp = W* uXp n^ l n x . 

Now, recalling that the g-unit normal to S is given by n M = Cl^ 1 ^ one obtains 

wp V = vr 1 ^, w* pv = sr 1 ?^. (30) 

From the definition (|27p one readily obtains 

Ppu — ^ Ppv: 

Wp uX p = (P W*)*^ A p - (P W)pv\p 

= ^{Pu[xw p ]p - Pp[xw p]v - n[xw p]s e s p U - n [p w* v]s e s \ p ). 

5.5.3 Regularity at X' 

In what follows, it will be shown that the components of the Weyl tensor W pu xp with respect to 
the frame v a given by 

= {Wpuxp o $K„« V V p s 

do not diverge as p — ¥ 0. In order to do this, one needs to discuss the expansions of n p , „a, 
Ppv, «V and 

We notice that the hypersurface to be considered is given in our coordinates by 

S = {4>{x^) = t - constant = 0}. 
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It follows that the normal vector and covector are given by 

Vfl 



pil-VtfPpcp") 1 / 3 ' 

vn 



n a = 0, 



P VQf3 a 



(1 - V 2 fl 2 f3 c f3 c ) 1 / 2 ' 



whence 



n t =p + 2Mp 2 + — p 3 (7M 4 + MM ab e a e b + (S a e a ) 2 ) + 0(p 4 ), 



n a = 0, 
h 
n 



* = p- 1 - 2M + ^JP(5M 4 - MM ab e a e b - (S a e a ) 2 ) + 0(p A ), 



-2p 3 e abc e b S c + 0(p 4 ). 

To compute v \ one notes the relations v \ = g p,p Z pv \ and that t pu \ — n a i a - f j_^\, where 

fcr^A = I det (g f j, v )\^r) rTl j, v x = -^il + 2Mp + 0(p 2 ))r} afJ , u \, 

with rja^^x the totally antisymmetric tensor of rank 4. Next one evaluates p M „ = g M „ — 2n^n v to 
obtain 

( v 2 n 2 (i + v 2 n 2 p c f3 c ) i_ 
p 2 (i-v 2 n 2 /3 c pc) p 2 

1 --2r>2r 1 



^v 2 n 2 p a 



\ 



-sir-lab + V*&P a P h ) 



from which 



Ptt = -P 2 - 4Mp 3 - — P *(13M* + MM ab e a e b + (S a e a ) 2 ) + 0(p 5 ) 
„2„ „b cc i r\l „3 



2^e a6c e 6 S c + O( (0 d ), 
/0" 2 ^a6 - — ^ (Af*(*a4 ~ e a e fc ) + 2M(M ab + S ab M cd e c e d - 2e {a M b)c e c ) 



~2(S a S b + 5 ab (S c e c ) 2 - 2e( a S b )S, 



O(p). 



In order to calculate expansions for the electric and magnetic parts and one makes 
use of the expressions (|3T)f and (|2T)]) . From the formulae for conformal transformations one has 



= rab + n~ l D a D b n + n~ 1 D c D c nh ab - 2n- 2 D c rtD c Qh ab 



and 



Xab = Q(Xab + ^h a b), 



(31) 



where £ = n M <9 M f2. Hence one obtains 



w ab = -Q 1 (r b + fi ^oDftfi + ft 1 D c D c nh ab -Xc C Xab + XcaXb c 

+n- 1 t Xab + ci~ l t Xc c Kb - 2n- 2 D c nDMab - 2fi- 2 t 2 h ab ) . (32) 



The seemingly more singular terms in the last expression can be eliminated using the Hamiltonian 
constraint 

= r-(x c c ) 2 +Xc d X cd 

= n 2 f + 4nD c D c n - qd c vld c vl - n 2 (x c c ) 2 + n 2 Xc dx cd + ^^x c c - 6E 2 , 
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so that equation (|32j) transforms into 

w ab = -jr 1 ff ab - \fh ab + jr 1 (D a D b n - \D c DM ab ) 

-fc c - n-^ixab - \xc c K b ) + xcaxt c - \xcdx cd K b ) ■ (33) 

To calculate w* b one needs to use the conformal transformation for the derivative of a (0, 2)-tensor: 

DaXbc = DaXbc + ^(Zda&Xbc + d^Xac + O^Xha ~ h ab h de d e Q.Xdc ~ h ac h de 8 'd^Xdb) ■ 

The latter, together with e abc = " fabc finally yield 

w* ab = -Cl- l D cXd{a e b) cd . (34) 

In order to compute and w*^ from the expressions (f3"3")) and one uses that n^c^ = 
and ri^c* = so that n^w^ — and n^w* — 0. With this and the symmetry of the tensors 
one obtains 

n a n b w ab _ n b w ba 
w t t = ,_ n2 , wta = r7— , 



— n, — h — * — h — * 

n a n°w* b n°w ba 



(n*) 2 ' 



w fn = -- 



In order to evaluate the formula (l33l one makes use of 



Tab — d c T a c b — d a T c c b + T a c b T c d d — T c d b T d 



— (S ab - e a e b ) - (M A {5 ab + e a e b ) + M [M ab + 2S ab M c c - 3S ab M cd e c e a 



p2 ~u~», m 

+4e (a M b)c e c + S a S b + 2S ab S c S c - 35 ab (S c e c ) 2 + 4e [a S b) S c e c ) + 0{p), 



and 



D a D b n = d a d b n - r a c b d c n 

= ~e a eb + 4Me a e b + -^pp{^e{ a S b )S c e c + e a e b (S c e c ) 2 + M A (8 ab + ^e a e b ) 

+M(4M ab + e a e b M cd e c e d + Ae [a M b)c e c ) + 4S a S b ) + 0{p 2 ). 

To complete the analysis, one also needs expansions for Xab- To this end we make use of the 
tensor Xab and the results in |18) . More precisely, one has that 

u - 3 
Xab = Q 2 ?Pab, ^ab = -^e {a e b}cd S c e d + 0(p~ 2 ), 

so that 

Xab = 3pe {a e b)cd S c e d + 0(p 2 ). 

To get from Xab to X ab one uses the corresponding conformal transformation formulae to find 
that 

Xab = Ze {a e b ) cd S c e d + O(p). 

Finally, it is noticed that 

E = Ti^dpCl = 0(p 4 ), l a hc = P e a bc + 0(p 2 ) 

It follows then from formula (|3"3"|) that 

w ab = p~ 2 M(5 a b - 3e e 6 ) + 0(p°), 

w ta = 2Mp 2 e ab ce b S c + 0(p 3 ) = -Mp 2 (3 a + 0(p 3 ), 

w tt = 4Mp 6 (S a S a ~ (S a e a ) 2 ) + 0(p 7 ). 

25 



Similarly, one finds from (|34j) that 

Kb = lP~ l {($ab + e a e b )S c e c - ie (a S b) ) + 0(p°), 

wl = 3p 3 e abc e b S c S d e d + 0(p 4 ), 

w* tt = -6p 7 (S a S a - (S a e a ) 2 )S b e b + 0(p 8 ). 

In view of the previous discussion, one has all the ingredients to compute the leading terms 
of the expansions of the components W a p.ys ■ Due to the length of the calculation this has been 
done in a tensor manipulation program. One obtains the following: 

Lemma 4. The components, W a p-y8, of the Weyl tensor W pv \ p in the frame v a are regular (i.e. 
non-divergent) at I' . 

6 Stationary vacuum solutions near the cylinder at space- 
like infinity 

Once the regularity of the various conformal field at I' has been shown, the last step in our 
analysis is very similar to the discussion in |27j . The proof consists of several parts: first, one 
starts by giving explicitly a solution to the conformal geodesic equations on X'; in a second step 
a stability argument is used to show that the construction of the cylinder at spacelike infinity is 
regular in a neighbourhood of X; finally, one needs to show that the whole construction does not 
depend on the choice of conformal factor on the initial hypersurface. 

Remark. It is important to stress the differences in the regularity of static and stationary fields 
at spatial infinity. In the static case all relevant fields are analytic. By contrast, as shown in [18] . 
in the strictly stationary case the relevant fields are never analytic as functions of asymptotically 
Cartesian coordinates. However, as already seen, the stationary fields have an analytic expansion 
in terms of radial and angular coordinates. This is the type of coordinates used in both the 
general construction of the cylinder at spacelike infinity and in the extension discussed in Section 
[5j As a consequence, all the relevant fields are analytic in these coordinates. 

In what follows, the word analyticity will be used to describe analytic behaviour with respect 
to the radial and angular coordinates. 

6.1 Setting the conformal Gauss system 

We consider now the regular finite initial value problem at spacelike infinity for stationary data. 
For this, we make use of the initial hypersurface 

S = {t = 0}, 

and set the initial conditions on S that generate the desired conformal Gauss gauge system. 
6.1.1 Initial data for the canonical conformal factor 

The initial data for the conformal factor, 8*, is prescribed — cfr. equation ((HJ — by means of the 
function 

n = 2n\DMD a n\' 1/2 (35) 

so that 

e» = K^n. 

It follows that 

K = p + 0(p 2 ), 

e* = ^\b a nb a n\^ 2 = P + o(p 2 ). 

Hence, the conformal metric evaluated on S and the induced metric are given, respectively, by 

g^v = 0*5/1^1 h ab = Q 2 h ab — K~ 2 h ab . 
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6.1.2 Initial data for the tangent vector to the congruence of conformal geodesies 

Initial conditions for the tangent vector x = dx/dr, where r is the parameter of the conformal 
geodesic, are set to be 

x _L S, g(x,x) = l, on S. 

In order to implement the requirement of having x orthogonal to S we consider the stationary 
Killing vector £. One has that 

£ = — r-(«o + Vt) on 5, 
P 



where i>o and v± are the first two vectors of the frame v a discussed in Subsection 15.21 It can be 
readily checked that (£,d/5) = (£, dip A ) = on S so that £ _L S. In order to obtain the right 
normalisation one considers 

Thus, 

SI K 

x = jt(v + vi) = -(vq + «i) = X a v a , (36) 

e p 

where 

X a = 5 Q a + 5? + 0(p). (37) 
6.1.3 Initial data for the 1-form / 

The initial data for the 1-form / is chosen in agreement with condition ((SJ) so that 

(/, d T ) = 0, pull back of / to S = k dn. 

It follows then that (/, x) — 0. This property, together with the choice ((35)) for the function k 
— see equation ([TO) — give 

e = e*(i-T 2 ). (38) 

6.2 Relating the various conformal gauges 

The analysis performed in Section [S] was done in terms of the metric g^ v . Now, we proceed to 
discuss how this metric and its Levi-Civita connection, V, are related to the metric g^ v and its 
Levi-Civita connection, V. From the relation g M „ = 0, 2 g^ u it follows that 

g ^ = n 2 5 M „ = e%„, n = n^e. 

The relations between the physical Levi-Civita connection V, the unphysical Levi-Civita connec- 
tions V, V, and the Weyl connection V are given by 

V - V + S(f), 

V = V + S(/), 

V = V + S(f), 

V = V + 5(0 _1 de), 

V = V + 5(fi _1 dn), 

where 

/ = / + rrMn = / + e-^e - frMji. 

From the particular form of the conformal factor given by equation (|3"8]) one has that (d0, x) = 
on S. As SI is independent of t and d t is orthogonal S, it follows that (dO, x) = 0. Thus, (/, x) = 
on S. Finally, from 

n = ^ on S, (39) 

K 
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and observing that the pull-back of / to <S is given by p 1 dp, it follows that 

f = f a a a , with f a = -S a ° + S a 1 on S. (40) 

A property of conformal Gaussian systems is that (/, x) = in the whole of the spacetime. Hence, 
one has that 

n = n</,±> (41) 

along the conformal geodesies. This last equation, together with equation (13U1) . allows to deter- 
mine II if the contraction (/, x) is known. 

6.3 Solving the conformal geodesic equations 

In this section we provide a discussion of the conformal geodesic equations with respect to the 
metric g and of its solutions. A solution to these equations is given by a spacetime curve x(r) = 
(f (t), p(t), ^ a (t)) and a 1-form /(r) along the curve. If expressed in terms of the frame fields, 
v a , and coframe fields, a a , the functions involved in the conformal geodesic equations are the 
components g a /3, g a ^, ^a 1 p and L a p- These functions extend by analyticity through I' into 
a domain where p < 0. If one assumes such an extension, one obtains the so-called extended 
conformal geodesic equations. It is important to point out that the initial data on S are analytic. 
Therefore, one can consider these equations in a neighbourhood of I'. Moreover, it turns out 
that the restriction of the equations to 2° can be solved explicitly. The solution one obtains is 
universal in the sense that it is the same for all stationary solutions with non-vanishing mass. 
More precisely, one has the following lemma, whose proof is the same as that of Lemma 7.2 in 

Lemma 5. The solution to the restriction of the extended conformal geodesic equations to I' with 
the components L a p as given in Subsection \5.4\ and initial data x = (0, 0, ip A ), X a = So a + Si a 
and f a = —S a ° + 5a 1 is given by 

x(t) = (t(t),0^ A (t)) = (t,0^ A '), 

fo = ~T^—, h = l, fA=0. 
1 + T 

This solution extends by analyticity to a domain < r < 1 + 2e for some e > 0. The extension to 
X of the conformal factor II determined by equations (|41[) and (|39[) takes on I' the value H = 1 . 

The previous lemma not only gives precise information about the conformal geodesies ruling 
T' , but it also shows that these geodesies extend analytically beyond r = 1. These facts are, in 
turn, used to show that there exists a solution of the conformal geodesic equations near I', and 
that this solution extends for sufficiently large values of the parameter r. 

Consider a smooth extension of S into a range where p < such that (p,tp A ) extend as 
smooth coordinates. We denote this extension by S ex t- Now, if the extension S ex t\S is small 
enough, then the initial conditions for the conformal geodesic equations (1551) , (|37[) and (|4T)|) extend 
analytically to S ex t — the precise range of p is not required as long as it is small enough. Therefore, 
the conformal geodesic equations determine near S ex t an analytic congruence of solutions to the 
extended conformal geodesic equations. From Lemma [SJ and making use of well-known results of 
the theory of ordinary differential equations — see e.g. — it follows that, using the same e as 
in the lemma, there exists p# > such that for the initial data 

f = 0, p = p\ V>) = V^', with|p'|<p # , 



and what is implied at these points by equations (f36|) . (|37| and (|40l) . the solution of the extended 
conformal geodesic equations 

f = t(t,p / ^ a ), p = p(r >P '^ A ') ) i> A = ^ A (T,p'^ A '), f a = f a (T,p',Tp A '), (42) 

exist for the values 0<r<l + eof their natural parameter and the function LI is positive 
in the given range of p' and r. Moreover, the Jacobian of the map (r,p',ip A ) — > x m (t, p 1 ,tp A ) 
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takes the value 1 + f on I' and for sufficiently small p# > it does not vanish in the range 
< r < 1 + e, \p\ < p#. Therefore, the functions (T,p\ip A ) define a smooth coordinate system 
in a neighbourhood O' of I 1 in M! . In particular, the relation implies that the curves with 
p' > cross J? + for r = 1. The set O' contains the following special regions 

o'nJ r+ ' = {T = i,p'>o}, 
x' = {o < r < l,p' = o}, 

1+' = {r = l,p' = 0}, 

and it is ruled by conformal geodesies. As a consequence of this discussion, in terms of the 
frame fields v a and the coframe fields a a , the metric g, the connection coefficients, T a -y of the 
Weyl connection V and the components of the tensor fields L a p, f a , W a p-/S extend in the new 
coordinates as analytic fields to O' . 

Finally, the conformal geodesies on O' and the fields discussed in the previous paragraph can 
be used to implement the construction of the manifold Af as described in Section [3l This is done 
by solving linear ordinary differential equations along the conformal geodesies, with the given 
analytical initial data on S. One obtains the following lemma: 

Lemma 6. For stationary asymptotically flat initial data (as described in Section^ the construc- 
tion of Section [3] leads to a conformal representation of the stationary vacuum spacetime which, 
in a neighbourhood O C Af of the set I, is real analytic in the radial and angular coordinates. 



6.4 The conformal gauge for the initial data 

The conformal representation discussed in the previous lemma has made use of the 3-metric h 
and the conformal factor ft on S. It remains to be verified that the whole construction is robust 
with respect to rescalings on the conformal initial data of the form 

h h = $ 2 h, n^n = m, (43) 

where •& is an analytic, positive conformal factor. These rescalings correspond to a change in 
the conformal gauge, and imply a harmless change of the normal coordinates x a — > x al with 
x a '(0) — and an associated change e a — > e' a of the frame vector fields tangent to S, which 
will be propagated along the new conformal geodesies. It is necessary to understand how the 
congruence of conformal geodesies corresponding to Cl relates to the congruence corresponding to 

n. 

The conformal rescaling given by (|4"3"f also implies the transitions 

k — >• k~ — , 6* — > 6* = <^8*, 

where 



The function <j extends to S as an analytic function of (p,ip A ). From the initial conditions for 
the fi-congruence of conformal geodesies one gets that 

* = ? _1 i ) fs = fs+^dd-^ds, 

where the subscripts indicate the pull-back to S. It can be verified that 

x _L S, <d 2 g(x, x) = 1. 

In what follows, we consider the equations for the il-congruence in terms of g and its Levi-Civita 
connection V. As a result of the conformal invariance of conformal geodesies, it follows that the 
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spacetime curves do not change (as set points) writing the equations in this form. Furthermore, 
their parameter f remains unchanged. The 1-form is transformed according to 

/->/* = /-(ee)- 1 d(ee), 

and therefore 

</*,i)=0, /J = f § + iTMtf, on 5. 

If one expresses the 1-form /* in terms of the g-orthonormal frame e a satisfying eo -L S, one 
obtains 

/o=0, f* = f a + <dtf,e ), a = 1,2, 3. 

The fields x and /* are the initial data for the fi-congruence written in terms of g, e a and V. 
As <r — > 1 and (di9, e ) = O(p) as p — > 0, then 

->■ 9, i ->• i, /* ->■ /„ as p ->■ 0. (44) 

This means that the limits of the initial data for both congruences coincide on X ' , and therefore 
the corresponding curves are identical on X' . 

Now, one can go back to the arguments used to show the smoothness of the construction in 
terms of the f2-congruence and apply them to the fi-congruence. One concludes that in a certain 
neighbourhood O' C M! of X' the gauge and construction of Section [3] for the ^-congruence is as 
smooth an regular the one based on the il-congruence. This result is summarised in the following 
lemma: 

Lemma 7. For stationary asymptotically flat spacetimes the construction of the set X' is inde- 
pendent of the choice of conformal factor Q. The set X' coincides with the projection n'(X) of the 
cylinder at space-like infinity as defined in Section^ 

This concludes the proof of our Main Theorem — cfr. Subsection 11.61 

7 Conclusions 

The discussion in the previous section has shown that, for initial data sets which are stationary 
in the asymptotic region, the construction of the cylinder at spatial infinity is as regular as one 
would expect it to be. As a consequence, the solutions to the associated regular initial value 
problem at spatial infinity are regular at the critical sets X^ notwithstanding the degeneracy of 
a subset of the evolution equations at these sets. As the length of our analysis shows, this is by 
no means an obvious result, and it makes evident the delicate interplay between geometry and 
properties of differential equations that the conformal framework allows to resolve. Moreover, it 
brings to the forefront the special role played by stationary solutions in the class of solutions to 
the Einstein field equations admitting a smooth compactification at null infinity. 

It is worth pointing out that the analysis carried out in this article is essentially a space- 
time one. A proof of our main theorem that relies only on properties of stationary data and 
the conformal evolution would be, by necessity, much more complicated and would require an 
understanding of the structure of the conformal field equations that is not yet available. 

It is expected that our analysis will play an essential role in the construction of suitable 
non-time symmetric generalisations of the rigidity results for asymptotically simple spacetimes in 
[46l [45] . In this respect, it will also be of interest to obtain a parametrisation of a large class of 
initial data sets for which it is easy to recognise when the data is, in fact, stationary. This type 
of characterisations may well require the consideration of other properties of stationary solutions 
at spatial infinity which have not been touched upon here — most notably, whether stationary 
data sets satisfy some generalisation of the regularity conditions of [Ml 12] • 
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A Various expansions 

In this appendix we collect the expansions of some auxiliary quantities used in this article. 
Recalling that 0. = p^^-V^Vl one obtains 

Q = p + Mp 2 + J^ p 3(5M 4 + 2MM ab e a e b + 2(S a e a ) 2 ) + 0(p 4 ), 

d a Cl = e a + 2Me a p + J^p 2 (l5M 4 e a + 2Me a M bc e b e c 

+4MM ab e b + 2e a (S b e b ) 2 + 4S a (S b e b )) + 0(p 3 ), 
8 a d b n = - p (S ab - e a e b ) + 2MS ab + -L^p(fM 4 (6 ab + e a e b ) 

+M(2M ab + 4e (Q M b)c e c + (S ab - e a e b )M cd e c e d ) 

+2S a S b + 4e {a S b) S c e c + (S ab - e a e b ){S c e c ) 2 ) + 0(p 2 ). 

Also 

g ab d a Cld b n = -p 2 - 4Mp 3 - ^Lp 4 (9M 4 + 2MM ab e a e b + 2(S a e a ) 2 ) + 0(p 5 ). 
The Christoffel symbols of g are given by 

It follows that 

Tt a t = -\g ah d h gu 

= p s e a + 6Mp 4 e a + -?—p 5 (2QM 4 e a + e a MM bc e b e c 

+MM a b e b + e a (S b e b ) 2 + S a S b e b ) + 0(p 6 ), 
f t b a = \{g bt d a g tt + g bc d a g tc - g bc d c g ta ) 

= -p 3 (2S c e ca b + e c S d e cda e b + e c S d e cdh e a ) + 0(p 4 ), 
f ab = g ct d {a g b)t + \g cd {d a gbd + d b g ad - d d g ab ) 

= p- 1 ^ - 28 c {a e b) ) + ^-Lp( - M 4 (S ab e c + e a e b e c - 25 c {a e b) ) 
+2M(2e (a M b) c - e c M ab + 25 c {a M b)d e d - 2S ab M c d e d 
+5 ab e c M de e d e e - 2e (a e c M b)d e d ) + 2(2e [a S b) S c - e c S a S b 
+2S c {a S b) S d e d - 2S ab S c S d e d + S ab e c (S d e d ) 2 - 2e [a e c S b) S d e d )) + 0(p 2 ). 
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